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C^; It has been suggested that the interactions of energetic particles with the foamy structure 

of space-time thought to be generated by quantum-gravitational (QG) effects might violate 

r^, ' Lorentz invariance, so that they do not propagate at a universal speed of light. We consider 

^ I the limits that may be set on a linear or quadratic violation of Lorentz invariance in the 

^ \ propagation of energetic neutrinos, v/c= [1 ± {E / MyQci)] or [1 ± {E / MyQG2Y]-, using data 

lO . from supernova explosions and the OPERA long-baseline neutrino experiment. Using the 
CN ■ 

Q ■ SN1987a neutrino data from the Kamioka 11, 1MB and Baksan experiments, we set the limits 
^ I MyQGi > 2.7(2.5) X 10^° GeV for subluminal (superluminal) propagation, respectively, and 
00 ! M^QG2 > 4.6(4.1) X 10^ GeV at the 95 % confidence level. A future galactic supernova at 
^ , a distance of 10 kpc would have sensitivity to M^qgi > 2(4) x 10^^ GeV for subluminal 
' (superluminal) propagation, respectively, and MyQG2 > 2(4) x 10^ GeV. With the current 

^ ■ CNGS extraction spill length of 10.5 fis and with standard clock synchronization techniques, 
■ " " ' the sensitivity of the OPERA experiment would reach Mi^qgi ~ 7 x 10^ GeV (M^qg'2 ~ 
8 X 10^ GeV) after 5 years of nominal running. If the time structure of the SPS RF bunches 
within the extracted CNGS spills could be exploited, these figures would be significantly 
improved to M^qgi ~ 5 x 10^ GeV {M^qg2 ~ 4 x 10^ GeV). These results can be improved 
further if similar time resolution can be achieved with neutrino events occurring in the rock 
upstream of the OPERA detector: we find potential sensitivities to Mi^qgi ~ 4 x 10^ GeV 
and M^QG2 ~ 7 x 10^ GeV. 
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1 Introduction 



Neutrinos from astrophysical sources and long-baseline experiments are powerful probes of 
potential new physics. They have already been used to discover and measure the novel 
phenomena of neutrino oscillations, thereby establishing that neutrino have masses [1,2]. 
Long-baseline neutrino experiments have also been used to set limits on quantum dccoherence 
effects that might be induced by foamy fluctuations in the space-time background in some 
models of quantum gravity (QG) [3-6]. It has also been suggested that the space-time foam 
due to QG fluctuations might cause energetic particles to propagate at speeds different from 
the velocity of light, which would be approached only by low-energy massless particles [7,8]. 
Any deviation from the velocity of light at high energies might be either linear or quadratic, 
5v/c = {E/Mqgi) or {E/Mqg2)^, and might be either subluminal or superluminal. Such 
effects are, in principle, easily distinguishable from the effects of neutrino masses, since they 
depend differently on the energy E. 

There have been many tests of such Lorentz- violating effects on photon propagation from 
distant astrophysical objects such as gamma-ray bursters [9], pulsars [10] and active galactic 
nuclei [11]. These tests have looked for delays in the arrival times of energetic photons 
relative to low-energy photons, and their sensitivities improve with the distance of the source, 
the energies of the photons, the accuracy with which the arrival times of photons can be 
measured, and the fineness of the time structure of emissions at the astrophysical source. 
The sensitivities of these tests has reached M^qgi ^ 2x 10^'^ GeV and M^qg2 ~ 4 x 10^° GeV 
for linear and quadratic violations of Lorentz invariance, respectively. 

At least one QG model of space-time foam [12,13] suggests that Lorentz violation should 
be present only for particles without conserved internal quantum numbers, such as photons, 
and should be absent for particles with electric charges, such as electrons. Indeed, astro- 
physical data have been used to set very stringent limits on any Lorentz violation in electron 
propagation. However, these arguments do not apply to neutrinos, since they are known to 
oscillate, implying that lepton flavour quantum numbers arc not conserved. Moreover, neu- 
trinos are often thought to be Majorana particles, implying that the overall lepton number 
is also not conserved, in which case QG effects might also be present in neutrino propaga- 
tion [14]. It therefore becomes interesting to study experimentally the possibility of Lorentz 
violation in neutrino propagation. 

Experimental probes of Lorentz violation in neutrino propagation are hindered by the 
relative paucity of neutrino data from distant astrophysical sources, and require the ob- 
servation of narrow time structures in neutrino emissions. However, there has been one 
pioneering experimental study of possible Lorentz violation using the long-baseline MINOS 
experiment exposed to the NuMI neutrino beam from Fermilab, which found a range of 
neutrino velocities -2.4 x 10"^ < {v - c)/c < 12.6 x 10"^ allowed at the 99% C.L. [15]. 
Assuming an average neutrino energy of 3 GeV, and allowing for either linear or quadratic 
Lorentz violation: v/c=[l± {E/M^qgi)] or [1 ± {E/M^qg2)\ the MINOS result [15] cor- 
responds in the case of hnear Lorentz violation to Mj^qgi > 1(4) x 10^ GeV in the case of 
subluminal (superluminal) propagation, and in the case of quadratic Lorentz violation to 
M^QG2 > 600(250) GeV. 
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In this paper we first establish hmits on Lorentz violation using neutrino data from 
supernova 1987a, using data from the Kamioka II (KII) [16], Irvine-Michigan-Brookhaven 
(1MB) [17] and Baksan detectors [18]. We find M^qgi > 2.7(2.5) x 10^° GeV for subluminal 
(superluminal) propagation, respectively, and M^qg2 > 4.6(4.1) x 10^ GeV at the 95 % 
confidence level. These limits are already much more stringent than those established using 
the MINOS detector. We then assess the improved sensitivity to Lorentz violation that 
could be obtained if a galactic supernova at a distance of 10 kpc is observed using the Super- 
Kamiokande detector, estimating sensitivities to Mj^q^i > 2(4) x 10^^ GeV for subluminal 
(superluminal) propagation, respectively, and MyQG2 > 2(4) x 10^ GeV. All these results 
are obtained taking neutrino oscillation effects into account, and assuming that any Lorentz 
violation is flavour-independent 0. 

We also discuss the sensitivity to Lorentz violation of the OPERA experiment at the 
CNGS neutrino beam from CERN |. We recall that the CNGS beam cycle provides two 
fast-extracted proton spills lasting 10.5 /is each and separated by 50 ms, each containing 
2100 bunches with standard deviation 0.25 ns, separated from each other by the CERN SPS 
RF bucket structure of 5 ns [20]. The OPERA data-acquisition (DAQ) system is organized 
in such a way that each subdetector provides its data with a distributed time-stamp with 
a granularity of 10 ns. If a time-synchronization method conceptually similar to that of 
MINOS between the CERN neutrino extraction-magnet signal and the OPERA time-stamp 
were implemented, the sensitivity would be greater than that of MINOS. This because, even 
though the baseline between the source and the detector are the same and the spill lengths 
are similar, the neutrinos in the CNGS beam typically have higher energies than those in 
the NuMI beam. Exploiting this feature, on the basis of an optimized analaysis we estimate 
that after 5 years of running sensitivities using OPERA could reach M^qgi ~ 7 x 10^ GeV 
{M^QG2 ~ 8 X 10^ GeV) for subluminal (superluminal) propagation. 

Further improvements in sensitivity would result if one could exploit the RF bucket 
structure of the spill. Assuming that the arrival time of the neutrinos would be correlated 
with the RF bunch structure with a timing accuracy of, say, 1 ns, the sensitivity to Lorentz 
violation could be improved to M,^qgi ~ 5 x 10'' GeV {MyQG2 ~ 4 x lO'' GeV) for the 
linear and quadratic cases, respectively. These results could be improved significantly if 
neutrino events occurring in the rock upstream from OPERA could be included in the 
analysis. In this case, the sensitivities would become M^qgi ~ 4 x 10^ GeV and M^qg2 ~ 
7 X 10^ GeV. In the case of quadratic Lorentz violation, this sensitivity is better than that 
obtained from supernova 1987a, and even improves on the sensitivity possible with a future 
galactic supernova. 

^This is a strong condition on any model of Lorentz violation, that is imposed by the success of conven- 
tional neutrino oscillation phenomenology, which implies that flavour-dependent dispersion effects can be 
neglected in the analyses of MINOS and OPERA data. Such effects could in principle appear in neutrinos 
from supernovae, but would not affect the results presented below, which are essentially independent of 
oscillation hypotheses. 

^For previous discussions of searches for Lorentz violation in neutrino data, see [14,19]. 
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2 Limits on Lorentz Violation from Supernovae 



In this Section we discus the supernova mechanism and the abihty to test Lorentz violation 
via the detection of neutrinos created in this process. We then analyze the data from the 
supernova SN1987a, the first supernova from which neutrinos have been detected, giving 
bounds at the 95% C.L.. Then we simulate a possible future galactic supernova and discuss 
the potential of the next generation of neutrino detectors, represented by Super-Kamiokande 
(SK), to improve this bound. 

2.1 Review of Neutrino Emissions from Supernovae 

The detection of neutrinos from SN1987a in the Large Magellanic Cloud (LMC) remains 
a landmark in neutrino physics and astrophysics. Although only a handful of neutrinos 
were detected by the Kamiokande-II (KlI) [16], Irvine-Michigan-Brookhaven (1MB) [17] and 
Baksan [18] detectors, they provided direct evidence of the mechanism by which a star 
collapses, and the role played by neutrinos in this mechanism [2]. The numbers and energies 
of the neutrinos observed were consistent with the expected supernova energy release of a 
few times 10^'^ ergs via neutrinos with typical energies of tens of MeV. A future galactic 
supernova is expected to generate up to tens of thousands of events in a water-Cerenkov 
detector such as SK, which will clarify further theories of the supernova mechanism and of 
particle physics [21]. 

Current simulations reveal several distinct stages of neutrino emission [22-24]. During the 
early stage with a typical timescale of a few milliseconds, huge numbers of Vf. are produced 
via pe — > n^e, known as the neutronization peak. Despite the huge numbers of neutrinos 
produced, these are difficult to detect water-Cerenkov detectors, because the neutrinos pro- 
duced in this process are detected via scattering on electrons and (in the case of the electron 
neutrino) via interactions with Oxygen nuclei. At the energies of interest the cross section 
for detection is dominated by the charged- current interaction ne"*", which detects anti- 

electron neutrinos. During the later stages of the supernova explosion, all flavours of neutri- 
nos and antineutrinos are produced with approximate Fermi-Dirac spectra, that are charac- 
terized by different average energies for different neutrino species: ^E^^ — (10 — 12) MeV, 
(E'pg) = (12 — 18) MeV and (J^v^ = (15 — 28) MeV (where denotes z/^, and their respec- 
tive antiparticles), with total emitted energy fractions s^^ = (10 — 30)%, Sp^ = (10 — 30)%, 
e,^ = (10-20)% [23,24]. 

The neutrinos produced in the supernova pass from densities close to nuclear density in 
the core through to the approximate vacuum of interstellar space, and the interactions with 
this matter dominate the neutrino oscillations. The neutrinos become maximally mixed at 
Mikheev-Smirnov-Wolfenstein (MSW) resonances and to first approximation the nature of 
the oscillations can be determined by the properties of these resonances. The resonance 
condition is A — Am^ cos 29, where A is the matter potential, Am^ is the difference in mass 
squared and 9 is the mixing angle. For a typical density profile and composition of the 
supernova medium, and typical neutrino energies ,the matter potential is positive (negative) 
for neutrinos (antineutrinos). Assuming just the three Standard Model neutrinos, there are 
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Table 1: The oscillation probabilities for the normal and inverse hierarchies, including the 
effect of the spectral split (SS), where the resulting Ve and fluxes are F^^ = pF^^ + (1 —p)F'^^ 
and Fp^ = pF^^ + (1 — P)Ep^ respectively. 

two possible MSW resonances, corresponding to the solar and atmospheric mass-squared 
splittings [25-27]. We know from the solar and KamLAND data that Am^i = — ml is 
positive, and therefore the corresponding MSW resonance is in the neutrino sector [28]. 

However, the sign of Am^j is undetermined and therefore the corresponding resonance 
could be in either the neutrino or the antineutrino sector, corresponding to the two possible 
mass hierarchies, the normal (inverted) for a positive (negative) the resonance 

there is a probability of transitions between the mass eigenstates, known as 'level crossing'. 
If the width of the resonance is large compared to the neutrino oscillation length at the 
resonance then the level crossing probability is small and the resonance is adiabatic. On 
the other hand, if the width of the resonance is small compared to the neutrino oscillation 
length scale, then transitions between the mass eigenstates occur and the resonance is said 
to be non-adiabatic. Combining current simulations of the supernova and the value of the 
solar mixing angle we can determine that the the solar resonance is adiabatic [1]. However, 
the current limit on ^13 is insufficient to determine whether the atmospheric resonance is 
adiabatic or not: simulations indicate that if sin^ 26'i3 > 10~^ the resonance is adiabatic 
and if sin^ 2^13 < 10~^ the resonance is non-adiabatic. The oscillation probabilities for both 
hierarchies are given in Table [H 

In addition to these effects, recent work has shown that neutrino self-interactions can 
induce large, non-MSW flavour oscillations [29]. These occur at large neutrino densities, 
just outside the neutrinosphere. For the normal hierarchy these effects have little effect on 
the flavour oscillations, but for the inverted hierarchy with non-zero ^13 significant fiavour 
changes can occur. These effects result in a 'spectral split', in which the and spectra 
are simply swapped above a critical energy, while the entire spectra of the and are 
swapped. For the case where the flavour transformations have occurred before the MSW 
resonances the flavour transformations can be thought of as changing the initial spectra, 
whereas in the case of shallow density proflles this becomes more complicated. 

We note in addition that, as the shock wave inside the supernova passes through the 
atmospheric resonance, it can change it from adiabatic to non-adiabatic, resulting in a time 
dependence in the signal that we do not consider in this paper [30]. 
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2.2 Analysis Techniques 



As previously discussed, it has been suggested that QG effects may lead to Lorentz- violating 
modifications in the propagation of energetic particles, and hence to dispersive effects, specif- 
ically a non-trivial refractive index. These dispersive properties of the vacuum would lead 
to an energy dependence in the arrival times of neutrinos. 

Even in the absence of any detailed, analytic understanding of time structure of a neutrino 
signal from a supernova, one can exploit the observation that, since the neutrino events have 
a range of energies, an energy dependence of the neutrino velocity would spread out the 
arrival times, compared to the signal if there were no dispersive properties of the vacuum. 
Any data set comprising both the time and energy of each neutrino event can be analyzed by 
inverting the dispersion that would be caused by any hypothesized QG effect. The preferred 
value of the energy-dependence parameter would minimize the duration (time spread) of the 
supernova neutrino signal. 

Assuming either a linear or a quadratic form of Lorentz violation: v/c = [1 ± [E / AI^qgi)] 
or [1 ± [E / MyQG2)'^], a lower limit on Mi^qgi and M^qg2 may be obtained by requiring that 
the emission peak not be broadened significantly. A non-zero value of M~q^^ or M~qq2 might 
be indicated if it reduced significantly the duration (time spread) of the neutrino signal. The 
duration (time spread) of the neutrino signal can be quantified using different estimators 
depending mostly on the amount of available statistics and time profile of the data set, if 
applicable El- In the following, we outline two estimators for analyzing neutrino signals, that 
we use first to quantify the limits obtainable from the SN 1987a neutrino data and then the 
sensitivities that would be provided by a possible future galactic supernova signal. 

2.2.1 Minimal Dispersion (MD) Method 

We assume that the data set consists of a list of neutrino events with measured energies E 
and arrival times t such as that in Table 2. In the first method, we consider event lists with a 
relatively low number of events, that do not allow a reasonable time profile to be extracted. 
In this case we consider the time dispersion of the data set, quantified by 



where t is the time of each detected event. We then apply an energy-dependent time shift 
At = TiE'-, where n = L/cMIqqi, varying M^qgi so as remove any assumed dispersive effects. 

The 'correct' value of the time shift r/ should always compress the arrival times of the 
neutrino events. Any other ('uncorrect') value of ti would spread in time the events relative 
to the 'correct' shift. Therefore, the dispersion ([1]) can be considered as an estimator to 
measure the degree of 'compression' of the neutrino events in time. In the following, we 
first apply this MD method in a warm-up exercise to the data from SN 1987a, and later we 
exhibit in subsection 12.41 the typical behaviour of this estimator versus ti for hypothetical 
data from a possible future galactic supernova. 

Statistically poor event lists, such as that for SN1987a, the only one currently available in supernova 
neutrino astronomy, do not allow the time profile to be classified, because time binning is impractical and 
one cannot apply nonparametric statistical tests to unbinned data. 
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Evaluating the dispersion ([T]) one obtains 



„2 / /'j. _ ZTii 



{{t - nE^ - {t - nE^)Y) (2) 

= {f) - {tf - 2r, {{tE^) - {t){E^)) + {{E'^) - {EY) . (3) 

Therefore, the dispersion of the 'de- refracted' time distribution is minimized by the param- 
eter r™", defined by 

We can use dl]) for any data set to estimate the scale M^qgi at which Lorentz violation is 
manifest. However, there are uncertainties in the energy and time measurements, as well as 
statistical uncertainties in the estimation of the observables calculated from any given data 
set, compared to their true values. We estimate the statistical uncertainties of an observable 

X clS 

Stat . li^^) - {^y 



where N is the number of events, and x = E,t or some combination of both. In order to 
estimate the uncertainties in r^^^"", we use a Monte Carlo simulation to repeat the calculation 
of T;*"*" including the energy and statistical uncertainties. We then make a Gaussian fit and 
use it to quote best-fit parameters and errors. 



2.2.2 Energy Cost Function (ECF) Method 

This is a different analysis technique that is mostly applicable to event lists that are sta- 
tistically rich. This means that one can combine the neutrino events into a time profile 
exhibiting pulse features that can be distinguished (parametrically or nonparametrically) 
from a uniform distribution at high confidence level. 

For the analysis we first choose the most active (transient) part of the signal (ti;t2); as 
defined using a Kolmogorov-Smirnov (KS) statistic. The KS statistic is calculated using the 
difference between the cumulative distribution function (CDF) of the unbinned data and 
that of a uniform distribution. The KS statistic is defined as the time that elapses between 
the minimum and maximum of this difference @. Having chosen this window, we scan over 
its whole support the time distribution of all events, shifted by At = tiE'', and sum the 
energies of events in the window. This procedure is repeated for many values of ti, chosen so 
that the shifts At match the precision of the arrival-time measurements, thus defining the 
'energy cost function' (ECF). The maximum of the ECF indicates the value of ti that best 
recovers the signal, in the sense of maximizing its power (amount of energy in a window of a 
given time width t2 — ti). This procedure is then repeated for many Monte-Carlo (MC) data 
samples generated by applying to the measured neutrino energies the estimated Gaussian 
errors. A typical ECF for one particular MC realization as well as the typical distribution 

The most active part of the signal can be also chosen by fitting the binned time profle, but the non- 
parametric way we use to extract a feature is less dependent on the time profile. In the case of a multipulse 
structure of the time profile, several windows may be analized separately. 
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of the positions of the maxima of the ECFs for many enegy MC reahzations are illustrated 
in subsection 12.41 (see Fig. [5] and Fig. [6]). 

We perform this procedure for different energy weightings E^, where n=0,l,2, summing 
up either the numbers of events, the energies or the squares of the energies in the time 
window selected, so to optimize the errors placed on the scale of Lorentz violation. 

2.3 Data Analysis 

For the analysis of SN1987a we use the uncertainties in Table 2, which were taken from [31]. 
In the case of a possible galactic supernova, we consider the Super-Kamiokande (SK) water 
Cerenkov detector, and we use the detector properties given in [30,32], where the energy 
uncertainties are modelled as o"^*^*^ = \/EqE, where Eq = 0.22 MeV. We note that the 
uncertainties in the time measurements are in general much less than the statistical and 
energy uncertainties, and we therefore neglect them in our analysis. 
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Table 2: The measured neutrino data from SN1987a, where we have omitted the events 
identified previously as background, and in each data set we define t = 0.0s for the first 
event. 

2.3.1 SN1987a 

Neutrinos from SN1987a were detected in three detectors, KII, 1MB and Baksan. The times 
and energies of the events are given in Table 2. The minimum dispersion was calculated 
1000 times for each data set to include the smearing from uncertainties. As an example. 
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Fig. [T] shows this smearing for the KII data set. From these distributions we can determine 
the best fit and the error, which are summarized in Table [31 We analyze similarly the data 
from 1MB and Baksan. As there is an uncertainty in the relative time measurements of each 
detector, we analyze each data set independently using the minimal dispersion method and 
then combine them to quote the final best fit and error, as shown in Table [3l 
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Figure 1: The distribution Tmin of 1000 Monte Carlo simulations of the KII data on neutrinos 
from SN1987a, including the smearing due to energy uncertainties. 
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Error 
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-0.00643648 


0.101162 


-0.304 
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Table 3: The best fits to the SN1987a neutrino data obtained using the minimal dispersion 
method. 

On the basis of this combined analysis. Fig. [2] shows the region which is excluded by the 
SN1987a data. Taking the distance to the supernova as L = (51.3 ± 1.2) kpc, the scale at 
which Lorentz violation may enter the neutrino sector is constrained to be 

M^QGi > 2.7 X 10^° GeV or M^qgi > 2.5 x 10^° GeV (6) 

at the 95% CL. for the linear subluminal and superluminal models respectively. The corre- 
sponding limits for the quadratic models are 

M^QG2 > 4.6 X 10^ GeV or M^qg2 > 4.1 x 10^ GeV (7) 

at the 95% C.L. for the subluminal and superluminal versions, respectively. 
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Figure 2: The regions of parameter space excluded by SN1987a, for subluminal (dashing) and 
superluminal (black) linear (left) and quadratic (right) models. 

2.4 A Possible Future Galactic Supernova 

800 n — • — — — ' — — — — — • — — — ' — — — • — — — • — ^ 




2 4 6 8 10 

Time (s) 



Figure 3: The time distribution of events predicted by our Monte Carlo simulation for the 
case of subluminal Lorentz violation at the mass scales M — lO^'^GeV and M — lO^^GeV. 

The detection of a galactic supernova would provide improved sensitivity to the scale at 
which Lorentz violation might enter the neutrino sector, due to an increase in the number 
of neutrinos which would be detected. The number of events would also increase because 
the current neutrino detectors are larger than those used to detect neutrinos from SN1987a. 
However, these effects would be partially offset because ti (x L and therefore the time-energy 
shift will be reduced if, as expected, the supernova takes place within the galactic disc at 
a distance ~ 10 kpc, compared to SN1987a in the LMC at a distance of ~ 51 kpc. The 
increase in the number of neutrinos which are expected to be detected are also because the 
next supernova is expected to be closer to the earth than SN1987a. For definiteness, we use 
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here a Monte Carlo simulation of the Super-Kamiokande (SK) neutrino detector, but note 
that other neutrino detectors could also probe this physics [33]. Simulations estimate that 
the number of events detected in SK from a supernova at 10 kpc would be of the order of 
10,000 [21]. In order to analyze at what scales Lorentz violation could be probed by the 
detection of galactic supernova neutrinos, we made Monte Carlo simulations with various 
levels of linear and quadratic Lorentz violation. We used the energy spectra of neutrinos 
from the Livermore simulation [22], which is shown in Figure [Hand the detector properties 
given in [32]. 




J I I r I I I I I I I I I I I I I 



Time [sec] 



Figure 4: The neutrino energy spectra from the Livermore simulation [22]. 

We show in Fig. [3] results from our Monte Carlo simulation including both charged- 
current and neutral-current events for linear sublminal Lorentz violation at the energy scales 
MyQGi = lO^^Gel^ and M^qgi = lO^^GeV, including oscillations corresponding to the nor- 
mal hierarchy and assuming that the atmospheric resonance is adiabatic. The signal has 
spread out and shifted in time, as we would expect. This time shift is unobservable because 
it is shifted relative to the signal in the absence of Lorentz violation, which in practice cannot 
be measured. We have applied the MD and the maximal ECF methods with various energy 
weightings to the Monte Carlo data with Mqgi = 10^° GeV in order to estimate the level of 
Lorentz violation. 
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Fig. [5] shows the ECF for one reahzation of the energy-smeared sample obtained applying 
MC to the measured neutrino energies with the Gaussian errors expected from SK. It exhibits 
a clear maximum, whose position may be estimated by fitting it with a Gaussian profile in 
the peak vicinity. Fig. [6] shows the results of such fits to the ECFs constructed for the 1000 
energy-smeared realizations. From this distribution we can derive the prefered value of r/. 




-0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 

z (s MeV"^) 



Figure 5: The ECF linearly weighted with energy from one realization of the simmulated 
time profile Fig. with neutrino energies smeared by MC applying to the expecteed energy 
resolution of SK, for the case of linear energy depending neutrino velocity. 

The results are sumarized in Table HI where we have defined rhi = M^qgi/ M^Jqqi., where 
M^Jqqi is the true scale of Lorentz violation and M^qgi is that deduced from the analysis 
method. Comparing these results, we find that the maximal ECF technique has greater 
sensitivity than the MD method, and that the linear energy weighting has the greatest 
sensitivity among the ECF analyses. We therefore use this in the following. 

We have performed simulations for both the normal and inverted mass hierarchies, with 
and without the spectral splits caused by neutrino self-interactions, for the extreme cases 
Ph = 0.0 and Ph = 1-0, and analysed them using the ECF method. The corresponding 
results are summarized in Table where we see that Lorentz violation can be probed with 
similar sensitivity for all mass hierarchies. 

The top three rows of Table [6] show the results of our analysis for the linear cases M^qgi = 
(10^°, 10^^, 10^^) GeV, using the minimal ECF method with no energy weighting, and making 
linear and quadratic fits. We see that data from a future galactic supernova could place strong 
95% C.L. limits on the range of M^qgi if it is lower than 10^^ GeV. In the limit of negligible 
Lorentz violation (Mj.qci > 10^^ GeV), we find the lower hmits M^qgi > 2.2 x 
and MyQGi > 4.2 x lO^^Gel^ at the 95% C.L. for subluminal and superluminal models, 
respectively. The bottom three rows of Table [6] show the corresponding results for the 
quadratic cases MyQG2 = (10^-^, 10^, 10^-^) GeV, again using the minimal ECF method with 
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Figure 6: The distribution of the postions of the maximums from fits of ECFs like in Fig. 
of 1000 realizations of the simmulated time profile Fig. with neutrino energy smeared by 
MC. 



Method 


95% C.L. 


Minimal Dispersion (MD) 
ECF 0th order 
ECF 1st order 
ECF 2nd order 


0.60 < mi < 2.37 
0.90 < mi < 1.29 
0.88 < mi < 1.26 
0.87 < mi < 1.27 



Table 4: The 95% C.L. ranges of rhi = M^qgi/M^qqi obtained using the different dispersion 
methods and various energy weights for a Monte Carlo simulation of a possible future galactic 
supernova for P = 0.0, assuming the normal mass hierarchy and M^qgi = 10^^ CeV. 



no energy weighting. We see that data from a future galactic supernova could place strong 
95% C.L. limits on the range of M,jqg2 if it is lower than 10^ GeV. In the case of large 
MuQG2, we find the lower limits M^qg2 > 2.3 x 10^ GeV and Mj,qg2 > 3.9 x 10^ GeV at the 
95% C.L. for subluminal and superluminal models, respectively, in the quadratic case. 

Although the ECF method is more sensitive than the MD method, it is not applicable 
to a statistically poor data set. The ECF method is best for the analysis of a feature in a 
distribution superposed on a uniform background, and the extraction procedure is possible 
only with a relatively representative (i.e., large) sample of events. This is demonstrated by 
simmulating a possible future extra galactic supernova which might take place at a distance 
similar to that of SN 1987a. The simulation has been performed in such a way as to have 
a sample with sufficient statistics to claim at least a 3-0" detection of Lorentz invariance in 
neutrino propagation. This would need about 600 events for the linear case, corresponding, 
assuming the sensitivity of SK, to a supernova at a distance of about 40 kpc from the Earth. 
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Figure 7: The time profile of 600 events, which could he detected by SK from a future extra 
galactic supernovae occured at the distance 40 kpc from the earth. The events are simmulated 
with respect to the energy spectrum given in Fig. with linear energy depending propagation 
effect, encoded at the level ti = 5.5 s ■ MeV~^ . 

An expected time profile is presented in Fig. [71 The signal Fig. [7] contains 600 events and 
the time distribution encodes a linearly energy-dependent propagation effect at the level of 
Ti = 5.5 s ■ MeV~^ , corresponding to M^qgi = 7 x 10^ GeV . This distribution does not 
demonstrate any significant feature that one could extract in a time window to be analyzed 
using the ECF. Therefore, we apply the MD method, which is better for a signal with poor 
statistics. The typical behaviour of the dispersion ([2]) versus r for one realization of the 
energy-smeared sample of the 600 simulated events is presented in Fig. [HI The distribution 
of Tmin given by (jl]) of 1000 MC simulations similar to Fig. [1] recovers, in this case, the 
encoded signal ri = 5.5 s ■ MeV~^ {M^qgi = 7 x 10^ GeV) at the 3-cr level. A similar 
simulation for the quadratic case would require about 400 events, which would correspond 
to a SN at a distance of about 50 kpc from the Earth for the SK efficiency. A 3-cr signal could 
be recovered if the dispersion effect was at the level T2 = 0.1 s ■ MeV"~^, which corresponds 
to M^QG2 = 7 X 10^ GeV. 

The minimal 3-a discovery statistics, which amounts to 600 (400) events for linear 
(quadratic) energy dependence, is defined for the MD method by the uncertainty in the 
denominator of (jlj), which reads ~ 5/-\/iV A/\/N) for either the simulated events or 
events from SN 1987a, where is the number of detected events. This means, that for 
the statistics of SN 1987a, a Lorentz-violating signal could be detected only at about the 
1-a C.L., corresponding to the bounds obtained in the previous subsection. In the case of 
limited statistics like SN 1987a, It is possible to estimate similar limits on Lorentz violation 
without the full MD machinery used in 12.3.11 However, such an estimate would implicitly 
assumes that the dispersion of the initial signal at the source is known. One could rely on 
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Figure 8: The dispersion ^ versus r from one realization of the simmulated time profile 
Fig. with neutrino energies smeared by MC applying to the expecteed energy resolution of 
SK, for the case of linear energy depending neutrino velocity. 

computer simulations of a supernova explosion [22], but this would introduce an element of 
model-dependent information into the analysis. The methods considered here do not assume 
any knowladge on the true profile (spread) of the neutrino signal at the source: instead, they 
remove any propagation effect that may be encoded in the time profile. 

3 CNGS and the OPERA Experiment 

In this Section we discuss the sensitivities to Lorentz violation in neutrino propagation that 
could be provided by the OPERA experiment in the CNGS neutrino beam. We first discuss 
the sensitivity to Lorentz violation that could be obtained using the spill structure alone, 
without taking into account its bunch substructure. In a second step, we consider how this 
bunch substructure could be exploited to improve the sensitivity, which could be possible 
if the timing resolution currently expected for the OPERA detector could be improved 
significantly. 

We first recall some of the details of the pioneering analysis of the neutrino velocity in 
a long-baseline neutrino beam that has been published by the MINOS collaboration using 
the NuMI beam [15]. This analysis compared the absolute timings of the detected neutrino 
events in the near and far detectors. The arrival times in the near detector provide a 
direct measurement of the neutrino intensity time-profile, consisting of either 5 or 6 batches 
separated by short gaps within a 9.7 yus long spill. The near and far clocks were synchronized 
absolutely by means of Global Positioning Satellite (GPS) receivers. The resulting systematic 
error of ±64 ns was dominated by uncertainties in the delays in the optical fibres that ran 
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Mass hierarchy 


95% C.L. 


NH P = 0.0 
NH P = 1.0 
IH P = 0.0 
IH SS P = 0.0 
IH SS P = 1.0 


0.90 < mi < 1.29 
0.90 < mi < 1.28 
0.91 < mi < 1.26 
0.90 < mi < 1.27 
0.91 < mi < 1.28 



Table 5: The 95% C.L. for rhi = MyQci/ M^qqi obtained using the ECF method for a Monte 
Carlo simulation of a possible future galactic supernova, for the normal (NH) and inverted 
hierarchies (IH), and including the effect of a spectral split (SS), where P is the level-crossing 
probability, and NH P = 1.0 is equivalent to IH P = 1.0. 



Model 


95% C.L. 


M^QGi = lO^OGeV 
M^QGi = WGeV 
M^QGi = lO^^GeV 
M,QG2 = lO^-^GeV 
M,QG2 = lO'GeV 
M,QG2 = lO^-^Gel^ 


0.90 < mi < 1.29 
0.64 < mi < 1.93 
0.22 < mi, 0.42 < mf^'^" 

0.93 < m2 < 1.23 
0.65 > m2, 2.3 < mfP*^^ 
0.19 > m2, 0.72 < mfP'^^ 



Table 6: The 95% C.L. limits on M^qgi and MyQG2 obtained using the KS statistic and 
the ECF method, for subluminal Lorentz violation with certain input choices of M^qgi (top 
three rows) and MyQG2 (bottom three rows). We give the 95% C.L. limits for subluminal 
(superluminal) propagation as rhi (nf^^^^ ); if a limit for rff^^^'' is not given then superluminal 
propagation has been ruled out at the 95% C.L.. 



between the surface antennae and the underground detectors. Including the jitter of the 
two GPS clocks, the total relative time uncertainty was a =150 ns. This analysis measured 
{v - c)/c = (5.1 ± 2.9) X 10-5 at the 68% C.L., or -2.4 x 10"^ < {v - c)/c < 12.6 x 10"^ 
at the 99% C.L., at an average neutrino energy of 3 GeV [15]. In the case of linear Lorentz 
violation, this would correspond approximately to M^qgi > 1.2(4.2) x 10^ GeV in the case 
of subluminal (superluminal) propagation. 

3.1 CNGS Beam Characteristics 

The energy spectrum of the calculated CNGS z/^ flux is reproduced in Fig. [9l Its average 
neutrino energy is ~ 17 GeV, significantly higher than that of the NuMI beam. Since the 
CNGS baseline is almost identical with that the NuMI beam, this gives some advantage to 
OPERA, assuming that it can attain similar or better timing properties. We also recall that 
the CNGS beam is produced by extracting the SPS beam during spills of length 10.5 fis 
(10500 ns). Within each spill, the beam is extracted in 2100 bunches separated by 5 ns. 
Each individual spill has a 4 — cr duration of 2 ns, corresponding to a Gaussian RMS width 
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of 0.25 ns [20]. 
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Figure 9: The expected CNGS neutrino beam energy spectrum [20]. 



3.2 Spill Analysis 

We introduce a 'slicing estimator', based on the fact that if some energy-dependent time 
delay is encoded into the time structure of the spill by propagation of the neutrinos before 
detection, one should observe a systematic increase in the overall time delay of events as their 
energies grow. Therefore, we propose cutting the energy spectrum of the neutrino beam into 
a number of energy slices, and searching for a systematic delay in the mean arrival times of 
the events belonging to different energy slices that increases with the average energy of the 
slice. 

In order to illustrate this idea, we perform a simple exercise simulating the sensitivity of 
the slicing estimator for a time delay depending linearly on the neutrino energy: At = tE, 
assuming ^ 2 x 10^ charged- current events, as are expected to be observed in the 1.8 kton 
OPERA detector over 5 years of exposure time to the CNGS beam. We envisage superposing 
all the CNGS spills with a relative timing error 6t. Since each spill has 2100 bunches, we 
expect about 10 events on average due to each set of superposed bunches. As a starting-point, 
before incorporating the relative timing error, the timing of each event has been smeared 
using a Gaussian distribution with standard deviation 0.25 ns, reflecting the bunch spread. 
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We display in Fig. [10] a sample of events in our simulation, before incorporating the relative 
timing error and any delay in propagation. The 5 ns internal time structure of the spill is 
clearly visible. 
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Figure 10: A superposition of the production times of neutrinos in CNGS spills reflects the 
hunch structure of the CNGS beam [20]. 

We now incorporate the uncertainty in the relative timing of the bunch extraction and 
the detection of an event in the detector. The overall uncertainty has three components: an 
uncertainty in the extraction time relative to a standard clock at CERN, an uncertainty in 
the relative timing of clocks at CERN and the LNGS provided by the GPS system, and the 
uncertainty in the detector timing relative to a standard clock in the LNGS. With the current 
beam instrumentation, implementation of GPS and detector resolution, it is expected that 
this will be similar to that achieved by MINOS in the NuMI beam, namely ~ 100 ns. Such 
a timing error renders essentially invisible the internal bunch structure of the CNGS spill, 
which looks indistinguishable from a uniform distribution generated with the same statistics, 
as shown in the upper panel of Fig. [TTl 

We next demonstrate in the lower panel of Fig. [TT] the effect of a time delay during 
neutrino propagation at the level of r; = 5 ns/GeV, as would occur if M^qqi = 4.8 x 10^ GeV. 
This would correspond to a total delay ~ 100 ns at the average energy of the CNGS neutrino 
beam. We see clearly its smearing effect at the beginning and end of the spill, due to the 
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Figure 11: The time structure of events in the CNGS beam, including a 100 ns timing 
uncertainty without (upper panel) Lorentz violation in neutrino propagation, and (lower 
panel) with a linearly energy- dependent time delay during neutrino propagation at the level 
of T = 5 ns/GeV. 

later arrivals of the more energetic neutrinos. Our 'slicing estimator' aims to quantify this 
effect. 

We smear the events with an energy resolution of 20%, and then cut the sample into 
slices of about 1000 events each with increasing energies. The asterisks in Fig. [12] show the 
mean arrival times of each slice, relative to the mean time of the superposed spills, using 
one particular smearing of the timing with a Gaussian error 5t = 100 ns. The triangles in 
Fig. [121 oil the other hand, show the mean arrival times of events in each energy slice if the 
propagation delays caused by an assumed value of r = 5 ns/GeV are included. We see clear 
differences between the asterisks and the red and triangles, that increase with the energies 
of the slices. 

By making many realizations of the event sample with the Gaussian 6t = 100 ns smearing, 
one can understand the significance of the shifts in the mean positions of the slices. Fig. [T3] 
shows the energy dependence of the shifts in the mean timings of the slices of 1000 events 
with a delay r/ = 5 ns/GeV encoded. These points may be fitted to a straight line 

A{t)=ri{E) + b. (8) 
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Figure 12: The mean arrival times of 1000-event slices with increasing energies without 
Lorentz violation in the neutrino propagation (asterisks) and with the effect of a time delay 
during neutrino propagation at the level of t = b ns/GeV (triangles) . The latter corresponds 
to M^QGi = 4.8 X 10^ GeV. One particular simulation of the OPERA experiment is shown: 
others are similar, exhibiting the expected statistical fluctuations. 

In general, when choosing the number of events for each shce, one has to strike a balance 
between the statistics of each subsample (which determines the precision of the determination 
of the mean arrival time of each slice), and the number of subsamples to be included in the 
fit. We choose the statistics of each slice so as to give comparable error bars for each energy 
bin. The propagation effect of interest to us is reflected in the slope n, while the intercept is 
an overall shift that has no physical significance. The sensitivity of the experiment to linear 
Lorentz violation at, say, the 95% confidence level (CL.) may be estimated by finding the 
value of the parameter ti which yields a fitted slope parameter that differs from a horizontal 
line (t; = 0) by 1.95a or more. We show in Fig. [TUthe confidence contours corresponding 
to 68% , 95% and 99% sensitivity levels in the (b,Ti) plane. From the upper (lower) edge 
of the corresponding ellipse one obtains ri95% = 4.9(2.6) ns/GeV at the 95% CL. for the 
subluminal (superluminal) propagation schemes, corresponding via 

M^QOi = ^rfi = 2.4 X 10^ [^^^-^f^) GeV (9) 
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to values of the linear Lorentz- violating scale M^ogi = 4.9(9.2) x 10^ GeV for the subluminal 
(superluminal) case, yielding a mean sensitivity □ to M^qgi ^ 7 x 10^ GeV. It is important 
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Figure 13: The measured shifts in the average arrival times of neutrinos in 1000-event slices 
with increasing energies, assuming a time delay during neutrino propagation at the level of 
r = 5 ns/GeV. 

to note that the slope and intercept are anticorrelated in such a fit, as shown in Fig. [TH Our 
conservative estimate of the limits corresponds to the upper (lower) edges of the elhpse. 

If the velocity of the neutrino depends quadratically on the energy of the neutrino, the 
slices should obey a parabolic fit 



A(t) =r,(^)2 + c. 



(10) 



Since the CNGS spill is in principle time symmetric, the estimated sensitivities for sub- and superluminal 
propagation should be the same. The difference between these numbers reflects the finite size of the simulated 
sample. Here and subsequently we quote the means of our sum and superluminal limits as estimates of the 
CNGS sensitivity. 
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Figure 14: The 68% (dashed dotted line), 95% (dashed line) and 99% (solid line) sensitivity 
contours for the case of linear energy -dependent fit 

Here the propagation effect of interest is parameterized by r^, wliile an overall shift is reflected 
in the constant c. The sensitivity contours at 68% , 95% and 99% CL are presented in Fig.fTSl 
According to the formula 

T^^ /-^CNGS \ a in3 ^-'^ GeV'^ 

M^QG2 = y — - — = 1.6 X lOM GeV, (11) 

after substituting rq95% = 0.066(0.022), we obtain M^qgi = 6.2(11) x 10^ GeV ^8x10=^ GeV. 

The stability of the slicing estimator has been checked by generating several data sets 
that have linear or quadratic dispersion effects artificially encoded. To test our level of 
sensitivity, we chose the Lorentz-violating parameters to be close to our estimations of the 
levels of sensitivities in the case of where dispersion effects are absent. The encoded values 
have been recovered for the linear ([8]) and quadratic ( fTOj) fits to slices containing the same 
numbers of events. Slight variations in the numbers of events in the individual slices do not 
change substantially the levels of sensitivity estimated for 1000-event bins. Another check 
has been performed using the minimal dispersion method described in Section 2.2.1. This 
method has been applied to the whole sample of about 2 x 10^ events expected to occur in 
the rock upstream of the OPERA detector, and results very similar to those of the slicing 
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Figure 15: The same as in Fig. calculated for the sensitivity of the quadratically energy- 
dependent fit ( fT^) . 



estimator have been obtained. Although the whole data sample is very rich statistically, the 
time distribution, given the 100 ns time uncertanty assumed in the current analysis, is still 
featureless apart the edges of the spill |^. 

Another check has been made analyzing the distortion of the shape of the spill at its edges. 
For this purpose, we analyze two independent histograms of the type shown in Fig. [TTJ One 
of the histograms is treated as a reference, while the other was shifted by introducing a 
time delay ri(q) for every event, corresponding to the linear (quadratic) propagation scheme. 
The shifted histogram has been compared to the reference (unshifted) histogram, and the 
parameter ri(q) increased until the difference between two histograms reaches the 95% CL. 
We find that this edge-fitting method has a factor 5 less sensitivity than that obtained earlier 
with the slicing estimator or the MD method. 

We recall that the OPERA detector may also be used to measure the arrival times of 
muons from 2 x 10^ neutrino events in the rock upstream of the detector. Information on the 
neutrino energy is missing in this measurement. Therefore, one cannot employ methods in- 
volving time-energy correlation information such as the slicing estimator. Methods requiring 
an energy-dependent time shift of the data, like the MD method, are also not applicable in 
this case, again because events in the rock do not have measured energies. Nevertheless, one 



^For this reason, the ECF technique described in Section 2 is inapphcable. 
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Figure 16: The left (right) spill edges fitted using 20000 detector events for scenario with a lin- 
ear energy dependence of the neutrino velocity. The solid red line is the reference histogram, 
while the points represent the shifted data. The solid black line represents the probability 
distribution function. 



can use methods that compare overall the time shift of the simulated data to the measured 
time distribution of the rock events. In this spirit, applying to the 2 x 10^ expected rock 
events the edge-fitting procedure described in the previous paragraph, we find a sensitivity to 
MyQGi ~ 2.4 X 10^ GeV, about three times better than previously, for the sensitivity to linear 
energy dependence, and the same level of sensitivity for the quadratic energy dependence. 

One can also modify the MD method for analyzing rock events. Namely, one could 
generate the reference spill and introduce an energy-dependent shift via the parameter ri(q) 
so as to make the dispersion of the shifted reference spill match as closely as possible the 
dispersion of the events measured in the rocks. However, due to statistical uncertainties 
the dispersion of each reference spill will be different to the dispersion of the rock events. 
If this uncertainty is much less than the increase in the dispersion of the rock events due 
to Lorentz violation then this method can be used. However, this limits the sensitivity to 
MyQGi ~ 3 X 10^ GeV for the linear propagation scheme, which is not as sensitive to other 
methods we have described above. From the other side the sensitivity of this modified MD 
method approachs to MyQG2 — 7 x 10^ GeV, which is similar to the slicing estimator. 

3.3 Bunch Analysis 

We now explore the additional sensitivity that OPERA could obtain if it could achieve a 
correlation between the SPS RF bunch structure and the detector at the nanosecond level. 
Sub-ns resolution could be obtained in OPERA with the help of additional specialized tim- 
ing detectors such as TOF hodoscopes 0. However, synchronizing the SPS and OPERA 

^We point out that it is sufficient to refer all measured far times to a well-defined plane perpendicular to 
the beam axis. 
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clocks with such a precision over a period of 5 years is a challenging task. With the new 
IEEE Standard Precision Time Protocol (PTP) IEEE1588 [34] it is possible achieve time 
synchronization in the range of 100 ns on an Ethernet network but not better; GPS clock 
synchronization at the ns level is also highly demanding. Standard 'One-Way' GPS tech- 
niques [35] can reach a precision of ~ 20 ns at best. Devices known as GPS disciplined 
oscillations (GPSDO) [36], containing high-quality temperature-controlled local oscillators, 
steered to agree with the onboard oscillators of the GPS satellites, can provide ultra-precise 
standard frequencies that could reproduce the CERN RF frequency. A more elegant but 
less standard method is called 'Common- View' GPS [35]: in this case two clocks (e.g., one 
at CERN and the other at LNGS) view simultaneously the same GPS satellite, thereby 
cancelling out the common errors (e.g., the satellite's local clock). It has been shown that 
the data recorded by the two GPS receivers can be processed offline to provide a timing 
uncertainty < 5 ns. Finally it has been shown that 'Carrier-Phase' GPS measurements [37], 
which use the carrier frequencies instead of the codes transmitted by the satellites, can 
achieve synchronization of clocks with uncertainties ~ 0.5 ns at the cost of extensive post 
processing. Turning to ground based solutions, the most precise atomic clock (the NIST-Fl 
used to define the UTC) has a long-term accuracy of 5 x 10^^^ or ~ 75 ns over 5 years. It 
would therefore not be sufflcient to bring two a priori synchronized clocks to the near and 
far locations to define the arrival times with the required long-term stability. Alternatively, 
next-generation accelerators, e.g., free electron lasers such as XFELs that aim to generate 
X-ray pulses with pulse durations down to tens of femtoseconds, will meet the challenge of 
finding new methods of ultra-stable timing stabilization, synchronization and distribution 
over several kilometres. These systems will most likely rely on optical timing synchroniza- 
tion. We can therefore imagine a phase-locked loop RF oscillator located at the far location 
remotely locked to the SPS RF system. These two systems would be connected and locked 
via stabilized optical fibre links U- To conclude, a combination of space- or ground-based 
solutions could probably provide the possible synchronization of the CNGS and OPERA 
clocks, and allow for systematic cross-checks to be performed. 

We now discuss how the sensitivity of the previous analysis could be improved by taking 
into account the 5-ns bunch structure of the CNGS spills. In Fig.[T7]we present one particular 
realization of a sample of simulated events which incorporates a relative timing error of Ins. 
Although the periodic bunch structure survives, the signal itself represents a time series 
with a relatively low signal-to-noise ratio. The latter implies that the proper deconvolution 
to extract isolated features cannot be made. In the other words, there is a problem in 
fitting the fine structure of the signal with an analytical function. Such a situation has 
been widely investigated and applied to the temporal profiles of gamma gamma ray bursters 
(GRBs) [38]. We therefore apply a cross correlation function (CCF) method similar to that 
described in [38] but differing only in details of its adaptation. Namely, we introduce the 

^We note that the temperature dependence of the refractive index of an optical fibre is typically 10^^/K, 
which corresponds to a drift of 5 ns for 1000 km and a temperature stability of 1°C. 
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Figure 17: A particular realization of the bunch structure with ^ 1 ns relative time uncer- 
tainty incorporated. The histogram is binned with a resolution sutable for resolving the bunch 
structure. 

temporal correlation of two time series A{t) and B{t -\- T;(g)) 



where A{t) is a Monte Carlo simulation of the events with no dispersion effects, B{t — tiE^) 
is the simulated data which has the time shift required to invert the effect of the energy- 
dependent dispersion, A{t) and B{t — tiE^) are the mean values of the corresponding time 
series, and (TA{t) and (TBit-nEi) are the standard deviations from these mean values. We 
average over several Monte Carlo simulations to include the statistical uncertainties as well 
as performing time and energy smearing due to the uncertainty in these measurements. 

We calculate the CCF(ri(q)) as a function of ti and find its maximum value. The value of 
T/ which maximizes the CCF is an estimate of the true value of r^. To find this estimate we 
fit a Gaussian to the peak of the resulting CCF function shown in Fig. [TSl Each realization 
produced an independent measurement of the CCF at a given value of the shift parameter. 
The process of iteration for every value of the shift parameter in Fig. [18] was repeated until 
the resulting distribution approached a normal distribution, which typically took about 100 
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Figure 18: The Gaussian fit to the CCF calculated for the case of a linear energy dependence 
with time smearing ~ 1 ns. 

runs. Using these normal distributions the values and the standard deviations (error bars) 
presented in Fig. [18] have been calculated. 

The sensitivity of the CCF can then be estimated by the precision of the position of 
the maximum for the Gaussian fit in Fig. [TSl For the case of linear energy dispersion, the 
maximum of the CCF is found at t{^^^ = —0.033 ± 0.036 ns/GeV if no time shift encoded 
in the simulated data. For superluminal propagation, when t\ > 0, one can estimate T[g5% = 
0.037 ns/GeV, which corresponds via ([9]) to M^qgi ~ 6.6 x 10^ GeV. For the subluminal 
case, one obtains T^gr,% = O.lns/GeV, which corresponds to M^qgi ~ 2.4 x 10'' GeV. The 
same CCF procedure may also be applied to the quadratic shown in Fig. [191 The 

limits deduced from the fit Fig. [I9] are M^qg2 = 3.6(4.9) x 10^ GeV ~ 4 x 10^ GeV. 

Repeating the CCF procedure for a time resolution above 2 ns, one observes no maximum 
correlation in a reasonable range of the shift parameter, as seen in Fig. [201 From this one 
can conclude that the bunch structure degenerates into an essentially uniform distribution as 
soon as the time resolution becomes bigger than ^ 2 ns, in which case the slicing estimator 
described in the previous subsection should be applied. 

If the same time resolution ~ 1 ns can be attained for events occurring in the rock 
upstream from the OPERA detector, the CCF method can also be used to analyze these 
data, which should amount to some 2 x 10^ events. We see in Fig. [2Tlthat the bunch structure 
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Figure 19: The same as in Fig. [721 /or the quadratic energy dependence. 

of the rock events is clearly visible if a time resolution 1 ns is achieved, despite the fact that 
the energies of the neutrinos colliding in the rock cannot be determined. The CCF calculated 
for the rock events is presented in Fig. [221 together with a Gaussian fit. The sensitivities to 
Lorentz violation now attain the levels of M^qgi = 4.3(3.2) x 10* GeV ^ 4 x 10* GeV for 
the linear case, and Mi,qg2 = 8.8(4.3) x 10^ GeV ^ 7 x 10^ GeV for the quadratic case. The 
sensitivity in the quadratic case is significantly better than the sensitivity estimated for a 
possible future galactic supernova. 



4 Conclusions 

We find from the SN1987a data lower limits on the scale of linear Lorentz violation in the 
neutrino sector, namely M^qgi > 2.68 x 10^° GeV and M^qgi > 2.51 x 10^° GeV at the 95% 
C.L. in the sub luminal and superluminal cases respectively. The corresponding limits for the 
quadratic models are M^qg2 > 4.62 x 10^ GeV and M^qg2 > 4.13 x 10^ GeV at the 95% C.L. 
in the subluminal and superluminal cases, respectively. We have also used a Monte Carlo 
simulation of a galactic supernova at 10 kpc to estimate how accurately Lorentz violation 
could be probed in the future. In such a case one would observe more events because 
of the larger fiducial volume of the SK detector compared to the previous generation of 
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Figure 20: The profile of the CCF calculated with a 2 ns time resolution for the case of linear 
energy dependence in neutrino propagation. 

detectors, and because the next observable supernova is likely to be inside the galaxy and 
hence closer than SN1987a. On the other hand, if the next supernova is closer than SN1987a 
then the energy- dependent time shift due to Lorentz violation will be reduced, reducing 
also the expected sensitivity. We performed simulations for both the normal and inverted 
mass hierarchies and for both an adiabatic and a non-adiabatic atmospheric resonance. In 
all scenarios it would be possible to probe Lorentz violation using the methods decribed in 
this paper. We used the minimal dispersion (MD) method and the maximal EOF method 
with a several energy weightings and have shown that using the latter with a linear energy 
weighting has the greatest sensitivity. Using this method we have shown that we could place 
limits up to M^QGi > 2.2 x 10^^ GeV and M^qgi > 4.2 x 10" GeV at the 95% C.L. for the 
subluminal and superluminal cases, respectively, for linear models of Lorentz violation, and 
M^QG2 > 2.3 X 10^ GeV and M^qg2 > 3.9 x 10^ GeV at the 95% C.L. for the subluminal 
and superluminal cases, respectively, for quadratic models of Lorentz violation. 

We have then explored the sensitivity to Lorentz violation in neutrino propagation that 
could be obtained using data from the OPERA detector in the CNGS beam. By comparison 
with the result already obtained by MINOS in the NuMI beam, OPERA would benefit from 
the higher energy of the CNGS beam, the larger statistics we assume, and the possibility 
of exploiting the bunch structure of the CNGS beam that we have explored. We find that. 
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Figure 21: A simulated realization of the bunch structure for rock events, incorporating a 
timing uncertainty ~ 1 ns. The histogram is binned with a resolution suitable for resolving 
the bunch structure. 

using standard clock synchronization techniques, the sensitivity of the OPERA experiment 
would reach M^qgi ~ 7 x 10^ GeV {M^qg2 

~ 8 X 10^ GeV) after 5 years of nominal running. 
If the time structure of the SPS RF bunches within the extracted CNGS spills of 10.5 fis 
could be exploited, which would require reducing the timing uncertainty to ~ 1 ns, these 
figures would be significantly improved to M^qgi ~ 5 x 10'' GeV {M^qg2 ~ 4 x 10"^ GeV). 
Using events in the rock upstream of OPERA, and again assuming a time resolution ~ 1 ns, 
the sensitivities to Lorentz violation could be further improved to M,yQGi ~ 4 x 10^ GeV for 
the linear case and MyQG2 =— 7 x 10^ GeV for the quadratic case. While still inferior to 
the sensitivity of the supernova limits in the linear case, the OPERA rock sensitivity in the 
quadratic case would exceed even that possible using data from a future galactic supernova. 
This and the fact that any accelerator limit benefits from better-understood experimental 
conditions would motivate the effort that would be required to achieve nanosecond time 
resolution for the OPERA/CNGS combination. 
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Figure 22: The CCF for rock events with time resolution ^ 1 ns in the case of linear energy 
dependence, compared with a Gaussian fit. 
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